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Newly developed high peak power lasers have opened the possibilities of driving coherent light 
sources operating with laser plasma accelerated beams and wave undulators. We speculate on the 
combination of these two concepts and show that the merging of the underlying technologies could 
lead to new and interesting possibilities to achieve truly compact, coherent radiator devices. 
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I. INTRODUCTION 

Free Electron Lasers (FELs) operate in many regions of the electromagnetic spectrum and have opened new and 
wider perspectives in the applied science; the technologies associated with ultra-small world and ultrafast phenomena 
are getting a great deal of benefit within such a context jl|. 

The large size and the non-negligible costs along with the fact that these devices allow few experimental stations, 
compared to third generation light source counterparts, the dedicated beam-time is less. The existing and planned 
facilities for short wavelengths exploit and foresee the use of large-scale accelerators and long undulators to generate 
light 0. An effort aimed at reducing both the size and the cost of FEL operating systems could provide a way of 
making these sources more accessible in smaller and more numerous laboratories as working tools. 

For this aim, a number of options have been explored. The efforts have been directed to the reduction of undulator 
length, using e.g. electromagnetic devices, hereafter referred as wave undulator (WU), or the size of accelerating 
system, by means of the Laser Plasma Acceleration (LPA) mechanism, and more compact undulators Q- 
In terms of the length, the first solution, that of the wave undulator (WU), replaces the standard magnetic undulator’s 
period that is on the scale of centimeters with that of an ’’electromagnetic undulator” on the micrometer scale. Using 
a WU, this reduces the length of the undulator by a factor of a thousand or more. The LPA option allows accelerating 
gradients of tens of GeV/rn instead of MeV/m, given by the conventional radio-frequency-based accelerating devices. 
This acceleration scheme allows a further reduction on the order of another thousand in terms of length. We must 
point out, however, that this solution demands for a laser driver that is not an insignificant device. 

The merging of the two solutions - a unique coherent photon generating device (undulator) and the LPA acceleration 
scheme - [4[ is, therefore, something worth exploring to drastically reduce the device dimensions and system cost. 
This paper contains an analysis of the operation of the specific case of FEL operating in the X-ray (energy and 
wavelength) region by using a LPA e-beam and a WU. With this analysis, we explore the potential of developing a 
high-gain FEL source scaled to more modest size. The application of both tools - LPA and WU - could be applied 
for many configurations of various e-beam energies and/or photon wavelengths. The variants we will add, however, 
with respect to previous proposals is the sheared beam geometry [&] and the use of the same laser system to provide 
both the accelerated e-beam and the electromagnetic undulator. 


II. COMPACT FEL: A FEW SUGGESTIONS FOR THE DESIGN CRITERIA 

We will proceed first by assessing the reliability of the PARSIFEL code jf|, which will be used in this paper to analyze 
the characteristics of the proposed device. This code, based on a combination of semi-analytical formulae and scaling 
relations, has been used in the past to analyze the behavior of different FEL facilities operating in various configurations 
such as oscillators 0, SASE (Self-Amplified Spontaneous Emission) [§), and nonlinear harmonic generation Q. It is 
routinely used as an on line bench-mark of the SPARC facility jlo}. Very recently an updated version, including the 


Corresponding Author; Electronic address: giuseppe.dattoli@enea.it 



2 


effects of the electron beam transport, has provided an important tool for the study of SASE FEL operation with a 
mismatched electron beam [To]. 

The comparison will be made by confronting the PARSIFEL predictions with previous results (either numerical or 
experimental), available in the literature. The code, albeit not comparable in terms of completeness with Genesis 
rill , Ginger pjJ, Prometeo [l3j, Perseo 0 or Medusa [3 has the advantage of being extremely fast and provides a 
clear understanding of the role of various physical parameters entering the laser process. For a full system design, of 
course, we would perform simulations with several codes as well in order to reduce the overall design risk. 

Once the reliability of the method has been proved, by obtaining already published results concerning the operation 
in either the oscillator or SASE, we will discuss the possibility of driving a SASE FEL with an e-beam provided by a 
LPA process injected into a WU. 

The idea of operating a SASE FEL, using a laser instead of an undulator, emerged since the very beginning of FEL 
Physics [Hj]- With the advent of powerful laser systems, it is no longer a conceptual proposal [13,13) an d significant 
efforts have been devoted towards such realization 0- Laser intensities in excess of AQ 2& W[m 2 have, indeed, been 
reported [3], the upcoming petawatt lasers aim at intensities of the order of 10 27 WJm 2 ^A\ and intensities of the 
order of 10 28 4 - 10 30 W/m 2 are foreseen at the Extreme-Light-Infrastructure (ELI) [3] and as described in the US 
Department of Energy’s Laser report [3] ■ Another candidate for operation with a WU is provided by the so called 
radio-frequency undulators (RFU) [3 , which have the potentials for going to periods around 1 mm, with reasonably 
large K values (comparable with, or even larger than, those achievable with a laser). 

The possibility of designing SASE WU FEL with simple scaling formulae has been recently discussed in [3|- We, 
therefore, remind that in this type of device, the undulator strength parameter, associated with a WU of intensity / 
and wavelength A, in practical units reads as: 
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and that the other quantities (Pierce parameter p , gain length L g and saturated power Pf ), characterizing the FEL 
SASE evolution, can be summarized as: 
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where the parameter p has a role of providing the signal growth rate and the transfer efficiency from the electron 
beam to FEL photons, L g is the gain length of the system and fixes the dimension (length) of the device. Pf is the 
final laser power at saturation. 

Finally, the SASE intensity growth vs the ’’undulator” length, in the longitudinal coordinate z , can be specified by 
the following logistic-like function [3 | : 
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In Fig[TJi, we report the power growth of a WU FEL operating at A = 1.35nm and using a CO 2 laser, as undulator 
driver, with a power corresponding to K = 0.3, (/ [^-] = 9 • 10 18 ). In Fig. [TJd, we compare the results obtained by 
the GENESIS and those from the eq[3j respectively. The comparison is good, at least for the present purposes. Any 
discrepancies are related to the differences in the ID and 3D representations of the interaction and also due to the 
fact that the GENESIS starts from the random noise in the electron beam. 

The previous analysis does not include, however, the longitudinal distribution of the optical pulse, which may be 
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FIG. 1: a) Power (W) versus WU longitudinal coordinate (m); b)Comparison between the prediction of GENESIS (triangles) 
and the semi-analytical procedure (boxes), each step corresponds to 5 ■ lCP 3 m. 


responsible for the corresponding dependence of the K parameter on the pulse profile. The problem should be 
carefully taken into account in the simulations. Just to give an idea of the importance of the problems involved in, 
we use a fairly simple argument useful to fix the operating parameters. 

If we assume that the optical pulse has a Gaussian shape with rms value sz we impose the following condition 


K 0 - K r 
K 0 
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(4) 


which states that over the pulse length, from the beginning of the FEL interaction to the maximum of the pulse, 
the inhomogeneous broadening effects, induced by the longitudinal intensity shapes, are within the gain bandwidth. 
From eq. 0] we find that 
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Imposing that Ti corresponds to 1/2 of the saturation length, we obtain the condition 
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( 6 ) 


If we accept a saturation length of 0.005m with a p = 5 • 10 -4 we should demand for a r.rn.s pulse duration oy « 15ps. 
If we require a K value around 0.3, the necessary intensity per pulse should bring a very high energy density (in 
excess of 10 8 J/m 2 ), which will eventually be realized in the next generation of high power lasers. 

LPA driven accelerators have made significant progress during the last years. Beams with energies comparable to 
those of conventional synchrotrons have been reached. Albeit the beam qualities are still insufficient, in terms of 
energy dispersion, furthermore shot-to-shot stability control is presently rather poor, it is however hoped that, in the 
near term, stable laser accelerated beams will be developped for applications. 

One of these, LPA based SASE FEL, has been proposed 


27], and its principle has been tested [28(. The preliminary 


results show that both radiation and tunability can be achieved from the beam emission inside the undulator. 

The set of parameters reported in Tab. U is not too far from those achievable in the near future. They have, 
therefore, been exploited in start-to-end simulations modeling LPA SASE FEL with the Genesis [29j]. In Fig|S] we 
have also included the PARSIFEL predictions, and the comparison yields more than a reasonable agreement. We can, 
therefore, consider the use of our semi-analytical approach, reliable for the purposes of a quick evaluation of the FEL 
performances within different configurations. 


E = electron beam energy= 400MeF 
Q = charge per bunch= 250 pC 
t =r.m.s. bunch length= 10/s 
e = Normalized Emittance = 1 mm ■ mrad 
\ u = undulator period= 1.5cm 
K =undulator strength parameter= 2 


TABLE I: Set of parameters used to simulate the LPA SASE FEL with Genesis. 
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FIG. 2: Power (W) versus longitudinal coordinate: comparison between the predictions of the Genesis (Triangles) and of the 
Analytical model (Boxes). Each step n corresponds approximately to 0.5m. 


To avoid misunderstandings and for the sake of completeness we add a further comparison (Fig. [3]) , relevant to the 
operation at 10 nm with the same beam parameters of Fig. [2J but with a slightly different undulator. The comparison 
is not that striking as in the previous cases, but still acceptable. 



FIG. 3: Same as Fig. [2j but for the operation at lOnm (each step is 5m) (E = 750 MeV, K = 1.9, \ u = 1.5 cm ). 


In the case of SASE FEL driven by LPA beams the electron beam is affected by a large e-beam energy spread, 
which, as it is well known, induces an inhomogeneous broadening causing a decrease in the gain and a consequent 
increase in the saturation length. Such a spoiling effect could be counteracted by a large value of the Pierce parameter 
that for the specific case reported in TablelH is p 0 — 5.85-10 -3 (the subscript o stands for homogeneous) the inclusion 
of the inhomogeneous broadening correction [261 ] 
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determines a reduction of more than 50% (p — 2 • 10 -3 ). For larger energy spread and lower beam current, different 
solutions should be devised to mitigate such a detrimental contribution. One possibility is the use of the so called 
transverse gradient undulators (TGU), which can both counteract the effect due to energy distribution and provide a 
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more efficient operation. In TGU devices a transverse tapering of the magnetic field is induced, providing an energy 
dispersion in the tapering direction (see FigHJ). 


X 



FIG. 4: Scheme of a transverse gradient undulator. 


To better clarify the working principle of the device we make reference to Fig[5]where we have reported the transverse 
profile of an undulator along with the associated tapering, realized by canting its upper and lower faces. 

The LPA produced e-beam are characterized by an energy transverse position correlation, therefore, the tapering 
should be realized in such a way that the portion of the beam with lower energy intercepts a region of the undulator 
with weaker field intensity and viceversa, so that all the bunch region radiate the same wave length, thus cancelling 
the inhomogeneous broadening effects. 

If the energy and the strength parameters are specified through the relations 


j(x) = 7 o (1 + -D 1 x) 
K(x) = (1 + kx) 


( 8 ) 


with D and n being the energy dispersion function and undulator gradient, respectively, we end up with the following 
condition necessary to cancel the inhomogeneous broadening effects 
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The presence of a non- zero dispersion function provides an increase of the transverse dimensions according to 
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Which is associated with a reduction of the energy spread (due to the conservation of transverse longitudinal phase 
space) according to 
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As we already remarked the energy spread induces an increase of the gain length according to the equation 
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FIG. 5: Scheme of a transverse gradient undulator. 


Therefore, by redefining the energy spread through eq. |TT]we should be able to quantify the beneficial effect of the 
transverse tapering. The beneficial effect of the reduction of the relative energy spread is counteracted by the increase 
of beam section, determining a decreasing of the current density, which in turn determines a reduction of the Pierce 
parameter as reported below 


P = Po 



(13) 


By including either effects in the definition of the gain length, we end up with 


L a 




- 1/3 

+ 0.64 




The gain length can accordingly be written as 
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In Fig. [6] we have reported L g /L° g vs. A, for different values of p e . 

The function exhibits a minimum, depending on the values of the ’’natural” inhomogeneous broadening, in correspon¬ 
dence of such a minimum, denoted by A* we can fix the optimum of the dispersion function as 


D* » —A* 

Ve 


(16) 


The minimum of the gain length occurs in correspondence of 

A* « v 70 - 333 ^ - 1 ( 17 ) 

which yields, for the dispersion function 

D* « —^0333^31 (18) 

for /i e 1 the eq. [18] reduce to 


D* » 1.15— 
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(19) 
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FIG. 6: Gain length vs. D for p e = 6(dot) and p e = 10 (continuous). 


The optimum gain length writes 


L 9 (A*,p £ )«l.l^L° (20) 

and can be exploited to get a first idea of how the various terms should be combined in the design of a FEL device 
operating with a large energy electron beam to mitigate the relevant negative effects. 

To illustrate the possible perverse effect of the energy spread on the gain length, we pick out, as reference parameters, 
those considered in ref. [30] and reported in TabHIl 


Undulator parameter 

K 

= 2 

Undulator period 

Au 

= 1cm 

Beam energy 

E 

= 1 GeV 

Resonant wavelength 

A f 

s lnm 

Peak current 

I p 

a 10 kA 

Energy spread 

cr e 

ps 10 -2 

Normalized emittances 


ps 1 pmrad 

Horizontal and vertical size 

U x 

, y ps 11.3pm 

FEL parameter 

P ? 

s 6 • 10” 3 


TABLE II: Plasma beam and FEL parameters. 


3l|). The line shift due to energy 
thus compensating the 


With the above list we find p £ ps 3.3, which determines an increase of the gain length with respect to the homogeneous 
case of about a factor 3. 

The use of the previous optimization criteria applied to the case of TabHIl determines an increase of the gain length 
of a factor 2 only, with an associated dispersion function D* ss 4 • 10” 3 m. 

The analysis developed so far is not dissimilar from that suggested in ref.[3t|, apart from unessential details due to 
the different parametrization formulae we used here (for further comments see ref. 
and undulator tapering are the same, but opposite in sign and cancel, therefore, each other 
effect of the energy spread [12] ■ 

A comparison between TGU and conventional undulator operation is shown in Fig. [7] where we have considered the 
same parameters of Fig. [2] 

Having sufficient reasons to be confident in our computational tool, we can discuss a scheme of a FEL wave undulator 
driven by a plasma accelerated e-beam. 

In the following, we will assume that the accelerating mechanism _relies on the so-called bubble accelerating regime. 
We use the scaling formulae derived by Pukhov and Gordienko J33| to fix the parameters of the LPA accelerating 
system. Accordingly, we will write the beam energy E e and the number of accelerated electrons N e in terms of the 






FIG. 7: FEL operation with and without TGU, the same parameters of Fig. [2] the transverse tapering parameters have been 
chosen to halve the effect of the inhomogeneous broadening (for further comments see the third of refs. 1321 ). 


laser parameters as: 


E P = 



N e 


kr e V P* 


( 21 ) 


where Pl is the laser power, P* = 8.5 GW is the threshold power for relativistic acceleration, k = 2 tt/X is the laser 
wave vector, <tl the length of the pulse and r e the electron classical radius. The parameter 5 depends on the energy 
distribution and on the energy spread, one can reasonably accept for a safe FEL operation. In the following, it will 
be treated as a free parameter. 

The acceleration length is linked to the laser parameters by the following equations: 
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Z r being the Rayleigh length and R the bubble radius, furthermore, the condition for the PLA requires cfl << R. 
In the following, we will assume that the electron bunch has the same length of the laser pulse, therefore, from the 
second of eqs. eh we find that the associated beam current is: 
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where Iq = 1.7 ■ 10 4 A is the Alfven current. The beam power can, therefore, be written as (P* = m ( c Iq): 

Pe — S ■ Pl 


(24) 


We can now use the previous formulae to get an idea of the numbers involved in the operation of a LPA WU SASE 
FEL. By keeping <j£ = yP, x < 1 we can link the radius of the accelerating bubble to the accelerating length according 
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to the identity: 
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Thus, finding for the laser power necessary to bring the system to the energy E e 
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Finally, the achievable peak current is: 
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We can now use the previous relations to derive the laser characteristics. We fix the following acceleration distance 
L acc = 2 • 10 4 A, (A = lpm), to reach a maximum beam energy of 200MeV, which, assuming X = 0.5, 5=1, yields 
the following set of parameters: 

2.6 • 10~ 5 m, P ^ 1.9 • 10 13 W, r L = 40/s, / = 1.03 • I 0 

The current is quite large, but such a value can be exploited to compensate the effect of the e-beam energy spread, 
which might be quite large. 

Let us now consider a WU provided by a laser with the characteristics used to derive the plots in Fig. [L] It is evident 
that if we use a head-on collision between such a laser and the PL produced beam, the scattered process determines 
photons in the range of gamma rays. 

If we use instead the sheared laser configuration suggested in ref. Tf and reported in Fig. [5] in this scattering 



FIG. 8: Sheared Configuration from ref. 04] 

geometry the electron and laser beams are nearly co-propagating, the induced Doppler shift is, therefore, reduced to 
yield soft, rather than hard X — rays or even 7 — rays. The scattered wavelength is, indeed: 

A(1 -p) ^ A F 

1 —/3cOs(0iab) sin(^Oiab) 2 
A 
4y 2 

If we use the parameters reported in Table Q] (A = 750nm and a shear angle with cos(9i a {,) = 0.9983), we obtain 
Ax — lAnm. The length of the interaction region can be fixed by the use of quite a simple argument. To have an 
estimation of the orders of magnitudes involved we note that the ’’equivalent undulator” period is: 

X SH ~ ^ 

u ~ 2 sin (%^) 


Ax = 
Af = 


(28) 
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(about 13/zm for the parameters reported in Table Hill) the corresponding gain length will, therefore, be expressed by 

\SH 

T _ _U 

9 ~ 47 rV3p 

by assuming that the saturation length is Ls = 20 L g and that p = 10 3 , we obtain an interaction length of less than 

0.012m. 

Regarding the Pierce parameter p we will write it in the form (valid for a circularly polarized mode) 


Po = 


K 2 ^\ 2 x r e niab 
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where ni a b is the e-beam density in the lab-frame, which is linked to the current by 

I 

Tllab — 7 > 

2j 7T CCT x@"y 

If we consider a round beam with 50 pm transverse length, we find for ni a b the value reported in Tab. m which, along 
with the other values therein reported, yields po = 7.14 • 10" 3 , and the inclusion of the effects of the inhomogeneous 
broadening yields a reduction of 50%. With these numbers we obtain the power growth reported in Fig. [9j in which 
the dot curve is relevant to the inclusion of other effects determining a further increase (another factor of 2) of the 
gain length. 


Al = laser wave length = 750 nm 

II = laser intensity = lO 22 ^- 

niab = lab-frame e-b density = 4.3 • 10 24 m -3 

K = 0.435 

e = normalized emittance = 1mm • mrad, 

\ u = equivalent undulator period = 13 pm 

TABLE III: 



FIG. 9: Power versus longitudinal coordinate (m) for a LPA WU FEL. 

An interesting aspect of this scheme is the possibility of exploiting one laser beam only, split into two parts, one 
being exploited to generate the PLA beam and the other acting as sheared WU. 

It is worth to note that, since K 2 oc I then we have p oc s/7, the SASE FEL power scales, therefore, as Pp oc oc 

I 4 /3 

The bare essence of the process we have described is, therefore, that of a parametric up conversion, in which the 
laser radiation is converted into that of shorter wavelength, by means of a non-linear device combining the PLA and 
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sheared laser beam undulator mechanisms. 

Before closing this paper we want to mention the possibility of exploiting a RFU undulator, namely a different flavor 
of WU. The present status allows the achievement of short periods (1.4cm) and large equivalent on axis magnetic 
field (on the order of 6 kG) {H}. They are powered by a microwave, provided by an X-band klystron, the field strength 
is associated with the injected power by the relation (35| : 

K 1 2 3 4 5 6 7 8 “ ^2 46 m] N- 2/3 P[MW] (29) 

The achievable values in terms of K values and of period length cannot, at moment, suggest application for fourth 
generation synchrotron radiation devices. They are, however, interesting for low gain oscillators and the experimental 
results of ref.[23j| support this possibility. One of the most appealing conclusions of [23| is the use of a configuration 
in which a beam of electrons undergoing a bunching in a first undulator is injected into a second undulator, where 
it undergoes a ’’seeded” coherent harmonic generation (SCHG). Even though the adjective ’’seeded” in not fully 
appropriated, since the emission is due to a harmonic bunching induced by the seeding in the first part of the 
undulator. The scheme is reminiscent of some theoretical proposals put forward in the past (36l. [jt! . The most 
significant breakthrough being, however, that the second undulator, where the emission actually occurs, is a RFU. 
The results reported in the experiment suggest that this undulator can be exploited to drive an oscillator (we will 
provide a detailed analysis elsewhere). We only note here that the achievable values of undulator strength and period 
may allow the operation in the UV region with an e-beam of modest energy. Such an oscillator could have interesting 
properties, associated with the possibility of piloting the RFU polarizations. The technical issues as the stabilization 
of K within a band of <C will be discussed in a separate publication. 

A further possibility is offered by these type of devices, if exploited together with other undulator devices. A promising 
scheme is the use of RFU undulator with low K value coupled with the second conventional undulator adjusted on 
an harmonics of the RFU. If a powerful seed laser, resonant with the RF undulator is exploited, the large seed power 
is capable of inducing bunching at higher harmonics, in spite of the low gain associated with the small K (< 0.1) 
value of the RFU. The harmonics supported by the conventional undulator can then grow due to the FEL mechanism. 
We have checked the feasibility of such a scheme, which is even more appealing because the RF undulator is not 
particularly demanding in terms of K-values and thus of power to be handled. A power of 6MW is, indeed, needed 
for RFU with K = 0.1, a wavelength of 10cm and N = 50. 


III. CONCLUSIONS 

In this paper, we have used simple analytical means based on (reliable) scaling formulae to analyze a number of 
possible solutions providing a kind of road map for the realization of compact short wavelength FEL source. From our 
analysis it emerges that pivoting elements of the discussion are plasma accelerated beams and laser wave undulators. 
On the other side, the idea of exploiting TGU undulator schemes is extremely interesting and should be pursued, even 
though a careful design modeling the effect of the sextupolar terms on the laser dynamics is needed. These terms 
may be responsible for an additional inhomogeneous broadening associated with the transverse electron distribution, 
which may nullify all the benefits induced by the TGU design. 

We believe that the use of RF undulators could provide an extremely helpful tool, also if exploited at the present state 
of the art. They may offer, indeed, a wealth of possibilities, including the superposition with an ordinary magnetic 
undulator, combined in such a way that a bi-harmonic geometry [38] is realized or to guarantee mild transverse and 
longitudinal tapering, as it will be shown elsewhere. 
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